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Abstract 

Relations between scattering and production amplitudes are studied in a microscopic 
multichannel model for meson-meson scattering, with coupling to confined quark-antiquark 
channels. Overlapping resonances and a proper threshold behaviour are treated exactly in 
the model. Under the spectator assumption, it is found that the two-particle production 
amplitude shares a common denominator with the elastic scattering amplitude, besides 
a numerator consisting of a linear combination of all elastic and some inelastic matrix 
elements. The coefficients in these linear combinations are shown to be generally complex. 
Finally, the standard operator expressions relating production and scattering amplitudes, 
viz. A = T/V and Q^m [A) = T*A, are fulfilled, while in the small-coupling limit the usual 
isobar model is recovered. 

1 Introduction 

In a very recent article [1] we have shown that several hadronic three-body decays of J /il), D and 
Dg mesons can be well described, up to moderately high energies, in a model for production pro- 
cesses derived from the so-called Resonance-Spectrum Expansion (RSE) [2]. The RSE formalism 
amounts to an effective description of non-exotic meson-meson scattering, based upon quark- 
antiquark pair creation and annihilation allowing transitions between an infinity of confined qq 
states and the meson-meson continuum. An essential feature of the RSE is that it gives rise 
to closed-form expressions for the S'-matrix and even the fully off-shell T-matrix. Hence, exact 
analyticity and unitarity properties, as well as a correct (sub) threshold behaviour, are manifestly 
satisfied. Moreover, the resulting meson-meson production amplitude can be derived exactly. 
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too, by summing the corresponding two-body Born series, the only assumption being that the 
third particle acts as a mere spectator [1] . 

In the present paper, we shall further develop the formalism introduced in Ref. [1], so as to 
cover the most general multichannel case in mesonic 3-body decays. Clearly, at higher energies 
competing inelastic 2-meson channels require that the production amplitude be described by 
a vector and not a scalar function. The underlying 2-body scattering T-matrix is then a true 
matrix. Furthermore, also the quark-antiquark sector needs an extension, as there can be mixing 
of different qq channels that couple to the same meson-meson channels. This is the case in e.g. 
the production of / = tttt and KK pairs, which both couple to the nn (= {uu + dd) / \/2) 
and ss channels, giving rise to the isoscalar scalar resonances /o(600) (alias a) and /o(980). 
Finally, having a general and exact — within the model assumptions — expression for production 
amplitudes at hand, we may carry out a detailed comparison with the ansatzes employed in other 
approaches, focusing on common features as well as clear differences. 

Under the spectator approximation, we assume that a pair of mesons is created out of one 
qq pair emerging from the original decay, accompanied by a non-interacting spectator meson. 
In this process, we only consider OZI-allowed [3] strong transitions to pairs of mesons. The 
intensities of transitions between the initial qq pair (its quantum numbers, including flavour, are 
abbreviated by a) and the various meson-meson pairs i allowed by quantum numbers are given 
by coupling constants gai determined via the recoupling scheme of Ref. [4]. In Sec. 3, we derive 
the matrix elements t{i —>■ u, E) of the scattering amplitude at total CM energy E = y/s, for 
transitions between the meson-meson channels i and u. In the same section, we establish a relation 
between the common denominator D{E) of all matrix elements ^(anything — ^ anything, E') and 
the numerators of diagonal matrix elements of t{E), the latter representing elastic scattering. 
The production amplitude a{a i,E), which is related to the probability of producing a meson- 
meson pair i, assuming that a qq pair emerges in the initial — here not described - stages of the 
decay process, is determined in Sec. 4. Note that the process giving rise to the initial qq pair 
plus the spectator meson can be either weak or strong yet OZI-suppressed (see Ref. [1] for some 
examples) . 

The central result of the present paper is a relation between the production and scattering 
amplitudes which can be formulated as 

a{a i,E) oc + i X"^ y2 {gai Xu{E) t{v iy,E) - g^u Xi{E) t{i v.E)} , (1) 

where the stand for momentum distributions that will be specified in Sec. 4 (Eq. 19). 

We thus obtain the result that the production amplitude is in the first place given by the 
common denominator of the scattering amplitudes. This implies that, within the RSE formalism, 
resonance poles are identical for production and scattering, at least in the spectator approxima- 
tion. Secondly, we find that the remainder of the production amplitude to the z-th two-meson 
channel is proportional to the sum of the differences between all possible elastic scattering am- 
plitudes t{y — > I/, i?) and the inelastic amplitudes for the i-th channel, t{i — *• v^E\ This does not 
spoil our conclusion about the resonance poles, since all T-matrix elements share the common 
denominator. Note, however, that the contribution of the term v = i vanishes in the expression 
between braces on the r.h.s. of Eq. (1). Consequently, the amplitude a{a i, E) for the produc- 
tion of a two-meson pair i does not carry any dependence on the amplitude for elastic scattering 
i ^ i. In those cases where the coupling to different meson-meson channels vanishes or can be 
neglected, implying a 1 x 1 T-matrix, the resulting production amplitude is solely determined by 
the common denominator D[E) [1]. 
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Some words are due about the ii'-matrix formalism. In general, and so also here, the T-matrix 
can be written as T = K/{\ — iK), where is a real symmetric matrix. So, at first sight, it 
seems that the common denominator of all T-matrix elements is given by 1 — iK. However, this 
is not the case. First of all, is a matrix and not just a real function. But even in the 1x1 case, 
where K can be represented by a real function, it has a denominator itself, the zeros of which 
are the /T-matrix poles, located at the real energies where some eigenphase shift passes through 
90°. The common denominator above is then the sum of the denominator of K plus —i times 
the numerator of K. In general, when i^' is a matrix, this relation involves the determinant of 
K. Hence, comparing i^'-matrix poles, lying on the real axis in the complex energy plane, and 
resonance poles, which are usually in the second Riemann sheet with respect to the nearest "open" 
threshold, is far from trivial. Moreover, for some resonances, like the a and the k (/^(^(SOO) [5]), 
the respective /T-matrix poles, corresponding to an elastic phase shift 6 = 90°, lie about 350-600 
MeV higher in energy that the real parts of the respective S'-matrix poles, while mixing with 
other resonances (/o(980) and Kq{1A30)) further complicates the picture. So rather than making 
ad hoc assumptions about poles of the production amplitude, we shall straightforwardly derive 
the latter, and then see what its properties are. 

A final remark here concerns Watson's [6] theorem for production. This theorem implic- 
itly relies on having a potential which is energy independent or only weakly energy dependent. 
However, this is not the case here, because the energy dependence of our effective meson-meson 
potential is far stronger than that of the scattering T-matrix. As a consequence, the energy 
dependence of the production experiment does not resemble at all the one of the T-matrix, and 
all exercises imposing the Watson "theorem" or theorems derived from it are inappropriate here. 
This issue is analysed in more detail in Ref. [7]. 

2 The Resonance- Spectrum Expansion (RSE) 

Scattering from a weakly coupled resonant source has been studied in a variety of different 
approaches. For such systems it is observed that resonances occur at energies that are close to 
the unperturbed spectrum of the resonant source. Widths and mass shifts can be determined by 
perturbative methods, and expressed in terms of pole positions of the resonances in the complex 
energy plane. 

Intuitively, however, perturbative methods do not offer the correct strategy for strong interac- 
tions. Since in the present paper we are interested in obtaining exact relations between scattering 
and production amplitudes, which are moreover based on a microscopic description in terms of 
quarks, we rather fall back upon an approximate yet exactly solvable theory or model. Such a 
manifestly unitary and analytic framework is provided by the RSE. 

The RSE aims at describing the scattering of meson pairs in non-exotic channels, thereby 
assuming that in the interaction region a meson pair may temporarely transform into a quark- 
antiquark pair through qq annihilation and subsequent creation. The transitions of the system, 
from meson-meson pairs to qq pairs and vice- versa, are described by an off-diagonal potential Vt 
in the RSE, linking these two sectors to each other. It has a maximum at an interquark distance 
To which depends on the average effective quark mass and runs from slightly less than 0.2 fm for 
bb to about 0.6 fm for light quarks. Furthermore, we assume that this mechanism gives rise to 
the dominant meson-meson interaction in non-exotic channels. Here, we limit ourselves to the 
case where Vt is considered the only interaction. 

The intermediate qq states are supposed to have an unperturbed confinement spectrum de- 
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pending on the quantum numbers of the system. Its energy eigenvalues are all contained in the 
two-meson scattering matrix, which, consequently, develops corresponding CDD resonance poles. 

An addtional nice feature of the RSE, which will turn out to be crucial for the construction 
of the production amplitude, is the possibility to obtain the closed-form scattering T-matrix 
both in configuration and in momentum space. In the former representation, coupled- channel 
Schrodinger equations with the proper boundary conditions directly lead to the full solution, 
while in the latter picture individual Born terms can be explicitly calculated and then exactly 
summed up owing to the general separable nature of the effective meson-meson potential, with no 
need to solve the Lippmann-Schwinger [8] integral equations. This allows to verify the correctness 
of the momentum-space approach in the scattering case, which is the only method at our disposal 
to describe production. As we shall see below, a similar Born series can then be written down 
and summed up. 

3 Scattering 

The building blocks of the RSE meson-meson scattering amplitude are the effective meson-meson 
potentials and the free two- meson propagators, which are graphically represented in Fig. 1. It is 
the philosophy of the RSE that confinement and decay can be separated. Hence, in the interaction 
region, a two-meson system can appear as a permanently confined system consisting of a valence 
quark and a valence antiquark. Possible intermediate crypto-exotic multiquark states are not 
considered in the RSE. 



V 



(effective 2-meson potential) (two-meson propagator) 

Figure 1: Graphical representation of the building blocks of the RSE two- meson scattering amplitude. 
The solid lines represent valence quarks and antiquarks as in usual Feynman diagrams. In contrast, 
the gray areas stand for all possible confining interactions, like gluon exchange, sea-quark loops and 
their higher orders. The effective meson-meson interaction is represented by V . Furthermore, although 
the mesons in the two-meson propagators are considered pointlike in the RSE, for clarity they are here 
represented by double quark lines connected by confining interactions. 

The dynamics of the intermediate qq states is described by a permanently confining Hamil- 
tonian He-, which has a complete set of eigenstates at eigenvalues representing the confinement 
spectrum. The other part of the strong interactions, generating transitions between a two-meson 
system and a qq state, is given by a transition potential Vt. Consequently, the effective meson- 
meson interaction in Fig. 1 is described by the operator 

V = Vt^ [E-H,f' Vt , (2) 

where E is the total invariant mass of the coupled-channel system. These interactions and the 
free two-meson propagators, both depicted in Fig. 1, can then be used in an ordinary Lippmann- 
Schwinger [8] approach to scattering. Nevertheless, in constructing the Born series, quark-loop 
contributions to all orders are automatically accounted for, as becomes clear from Fig. 2. 
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Figure 2: Graphical representation of the RSE scattering amphtude. 



The RSE scattering amphtude depicted in this figure has the Born term 
V{t ^ j) = (z, p^ I [E{p) - Vt \j, p, ) . 
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In Ref. [2] it was shown how, under the RSE assumptions, the integrations can be done analyt- 
ically. For the present discussion, it is only necessary to mention the generic form of the RSE 
expression, given by 



^2 oo 

+ 1) Pe {pi ■ p, ) jMro) JiiPjro) Z^f{E) 



(4) 



The overall coupling X and the interaction radius tq represent the total probability of quark-pair 
creation/annihilation and the average interquark distance at which such processes take place, 
respectively; je stands for the spherical Bessel function for the relative meson-meson angular 
momentum i; pi and pj are the relative linear momenta in the two-meson channels i and j, 
respectively. The matrix Z is a real and symmetric function of the total invariant mass of the 
system. 

The intermediate qq systems may have different orbital angular momenta, for the same quan- 
tum numbers. For example, a meson-meson system with J^^ = 1 couples to qq systems in 
either an S* or a D wave. On the other hand, isosinglet mesons can be mixtures of different 
quark flavours, usually nn and ss. In such cases, more than one type of qq states are involved in 
the quark loops of the process depicted in Fig. 2. However, since the recurrencies of the perma- 
nently conflned qq systems are automatically included by the deflnition of the Born term (4), the 
number of qq channels that couple to a speciflc set of two-meson quantum numbers is limited, 
usually to one or two. Nevertheless, the number of two-meson scattering channels is in principle 
not restricted. 

For any number of coupled conflnement and scattering channels, the general structure of the 
amplitude reads {E = E{pi) = E{pj)) 



t{z-.j) = {t,p, \t\j,p,) = {i,p, \{y + VGV + VGVGV + . . .)\3, p,) 
1) Pi {pi ■ Pj ) ji (piVo) je (pjTo) 



(5) 



^2 oo 

Air' 



A^{E) 



-I] 



where A and T) are functions of the total invariant mass E satisfying the unitarity condition 

2A2E/^^P^j1(P^^o)^l'U5f . (6) 



The denominator V contains the full pole structure of the coupled two-meson states. In order 
to be a bit more speciflc, let us consider the scattering of charmed mesons, i.e., DD, D*D, D*D*, 
DgDs, D*Ds and D^D*, all coupled to cc. For such a process, V has in the RSE the form 



= l + 2^^' 



n=0 



(ro) 

E — En 



l^uPuje (Puro) (p^ro) 



(7) 
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where the outer sum runs over all two-meson channels, and the inner sum over all recurrencies n 
for the operator He describing confinement in the cc system. F^"^ and En represent the eigenstate 
and eigenvalue of the n-th recurrency of the He spectrum, respectively. Furthermore, the g,^ stand 
for the relative couplings of each of the two-meson systems to cc, while hf^ is a spherical Hankel 
function of the first kind. 

The denominator T>{E) vanishes for E near En and small overall coupling A. In this case, 
the scattering cross sections in all channels display narrow spikes for values of E in the vicinity 
of En {n = 0, 1, 2, . . .). Hence, for small A, the theoretical cross sections reproduce — up to 
small shifts — the hypothetical cc confinement spectrum. 

However, for larger values of A the zeros in T) are no longer near the eigenvalues of He-, but 
move deeper into the complex E plane, farther away from the real axis and with appreciable shifts 
for the real parts as well. Then, the resonance spectrum does no longer reproduce the spectrum 
of He.: resonances start overlapping and even the number of zeros in T) that lie close enough to 
the real energy axis to be observed experimentally may change. We believe this describes quite 
accurately the true situation in hadron spectroscopy. 

Below the lowest threshold, poles, i.e., zeros in V (Eq. 7), come out on the real axis, because 
the expression ijihf^ turns real. In that case, expression (5) describes bound cc states, such 
as ?7c, J/il^i Xc(l-P) and ip{2S), yet with an admixture of two-meson components. The energy 
eigenvalues of these "dressed" states then depend on the value of A. It has been observed [9,10] 
that charmonium mass shifts with respect to the pure confinement spectrum can be surprisingly 
large in the RSE, as well as in other approaches [11]. 

In the present work, we intend to derive relations among A[f, V^^^ and zl-\ In principle, 
this could be achieved by just performing the calculus outlined in Ref. [2]. However, here we 
shall allow more general expressions for the Z matrix in the Born term (4). Hence, apart from 
the unitarity condition (6), we must construct a second relation. For that purpose, we write the 
identity 




+ 2 oo 

p,\{T-V- TGV) I J, p;. ) = — ^(2£ + 1) (p, ■ ) jeip.ro) n{p,ro) x 

which yields the relation 

4) = A? + 2^X'Y: l^Mi (Puro) hf (p.ro) z\^ . (9) 

Furthermore, if we assume 

A§ = Ai?^^ + AMlf^) + AMSf^) + ... , (10) 
then we obtain the following solution to relations (6) and (9): 

1. The denominator V can be fully expressed in terms of the numerators A, according to 

= l + 2zX'Y. f^'^P^Ji (P-^o) hP {p^ro) A^l . (11) 
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2. The zeroth- order term of (10) is evidently given by the Born term (4): 



J^"" = . (12) 



3. For the higher-order terms of the expansion (10) we obtain the recursion relation 

A^f = 2^Y. i^.V.U iP.ro) h'P (p.ro) {aIT' - Z^} . (13) 

V 



Prom Eq. (5) we then get a partial-wave scattering amplitude of the form 

te{i^j) = '2X'^je{piro)je{Pjro) :^W^ ■ (14) 

For a full definition of this amplitude, satisfying the unitarity conditions for scattering, see 
Eq. (29). 



4 Production 

Various opinions exist on how to analyse the final-state interactions of pairs of hadrons emerging 
from a decay process [12-15]. In particular, the production of pion pairs has been studied from 
many different angles. Several resonances have been discovered and established in this channel. 
However, there still are many open questions, of which the most intriguing one probably is the 
formation of the /o(980) resonance [16-28]. As such, this resonance seems to be one of the key 
issues for understanding strong interactions. It lies close the KK threshold, couples relatively 
weakly to pions, comes on top of a much broader structure, namely the /o(600), and is furthermore 
not very distant from a broad resonance around 1.35 GeV, viz. the /o(1370) [15]. 

It is our understanding that mesonic resonances, like the /o(600) and the /o(980), form 
an integral part of the whole meson family. Therefore, we have developed a model for all qq 
phenomena, including those involving charm and bottom. Here, we wish to develop a new tool 
for data analysis, which is an amplitude for the description of final-state interactions in two-meson 
subsystems emerging in decay processes involving other particles. This production amplitude is 
based on the two- meson scattering amplitude given in Eq. (5). 




+ ••• 



Figure 3: Graphical representation of the RSE production amplitude. The transition qq — > MM, 
denoted by Vt in the text, is here represented by v; the resulting effective MM interaction is denoted 
hyV. 

For the description of the final-state interactions of meson pairs in production processes, it 
is common practice to make the spectator assumption, according to which the other emerging 
hadrons do not interact strongly with the pair. Evidently, this is an approximation, which is 
justified by the observation that in most production processes resonances involving the third (or 
fourth, . . . ) hadron are much higher in mass than the energies considered for the pair. Here, 
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we moreover assume that the meson pair is generated from an initially produced qq pair. Our 
amplitude for the production of a meson pair, including all higher-order contributions from final- 
state interactions, is depicted in Fig. 3. Also using expression (5) for the scattering amplitude, 
we are led to define for the production amplitude 

a(« ^i) = {t, p1 1(1 + TG)Vt\ (qq)^ , E) = (15) 
= {h P^ \Vt\ (qq)^ .E) + Y.[ £K (t, Pi \T\ v, K)G (k) {i^, K \Vt\ (qq)^ , e) 



= ^ E i-^Yj^ fe^o) iP^ ) Q[:l (E) X 

Here, Q^"] represents the overlap with the initial qq distribution, having quantum numbers a and 
relative interquark angular momentum igg. Notice that the latter quantum number is related 
- though unequal - to the relative two-meson angular momentum i, because of total- angular- 
momentum and parity conservation. Below, we shall discuss the properties of production ampli- 
tude (15) for pairs of interacting mesons. 



4.1 Pi Cj/X'jy^? 

V 

The result (15) agrees to some extent with the expression proposed in Refs. [29,30]. Like here, 
the authors of Ref. [30] based their ansatz on the OZI rule [3] and the spectator picture, so 
as to find that the production amplitude can be written as a linear combination of the elastic 
ti{i — >■ %) and inelastic t^(i ^ v ^ i) scattering amplitudes, with coefficients that do not carry 
any singularities, but are rather supposed to depend smoothly on the total CM energy of the 
system. 

Indeed, if we carry out the substitution (14), we find for our production amplitude the ex- 
pression 

a{a ^i) = ^^(-z)^y"jf) {Pi)Qtl (E) IgaiJeiPiro) l^^Pu^f^ i.Pu'^o) dautS ^ , 

V ^ e,m I ) 

(16) 

which contains a linear combination of elements of the scattering amplitude, with coefficients 
smooth in E. 

However, Ref. [30] concluded from the relation 

^m{A) = T* A (17) 

that the production amplitude must be given by a real linear combination of the elements of the 
transition matrix. A similar conclusion, based on a i^-matrix parametrisation, can be found in 
Ref. [31]. In contrast, we arrive at a different conclusion, namely that, as the Hankel function 
of the first kind is a complex function for real arguments, the coefficients must be complex, in 
agreement with experimental analyses [13, 32, 33] as well as with the theoretical work of the 
Ishidas [34,35]. 
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Relation (17), which can be also found in Ref. [36] basically stems from the operator relations 
AV = (1 + TG)V = V + TGV = T, the symmetry of T, the realness of V and the unitarity 
of 1 + 2zT, which gives (A) V = (AV) = 3m (T) = T*T = T*AV. This leads, for non- 
singular potentials V, to relation (17). In Appendix A, we show that notwithstanding the complex 
coefficients in Eq. (16), relation (17) is satisfied for the scattering and production amplitudes of 
Eqs. (5) and (15), respectively. Consequently, relation (17) does not impose a realness condition 
on the coefficients in Eq. (16). 



4.2 The lowest-order term 

Besides the sum over transition matrix elements, our prodution amplitude (16) also contains 
an extra term oc Qai je (Pifo)- Such a term was not considered in Refs. [29-31]. However, in 
the works of Graves-Morris [37] and Aitchison & collaborators [38-40], the possible existence 
of an additional real contribution was anticipated. Here, it follows straightforwardly from the 
reasonable assumption that the produced meson pair originates from an initial qq pair. 

It is generally agreed that production and scattering have the same singularity structure in 
the complex energy plane. At first sight, this is not obvious from expressions (15) and (16). 
However, the second term between braces in Eq. (15) can, using Eq. (11), be rewritten as follows: 

g^i - 2iX^ l^uPi^ji (Pv^o) hf^ (p^ro) gau = (18) 



1 



9ai 



gai + 2iX^ fJ'uPujeiPuTo) (PiyTo) \gaiJ^yl - gauA 



■ iu 



2i Y l^yPvjf- {PvTo) /if ^ (Pi/ro) 



9o 



9o 



A 



((■) 



From this equation it is obvious that, in our approach, scattering and production have exactly 
the same poles in the complex energy plane, as they share the global denominator T). 



4.3 The central result 



The pole structure of our production amplitude is exhibited very explicitly in formula (18), and 
shows that it is completely given by P, the very same denominator that determines the pole struc- 
ture for elastic scattering. The conclusion is that resonance shapes are different for production 
and scattering because they are largely determined by the respective numerators. Moreover, pre- 
cisely the numerator An describing elastic scattering in the z-th two-meson channel has dropped 
out of expression (18). Hence, when restricted to a one-channel model, our production amplitude 
is completely determined by just the denominator V. 

The result (18) may be substituted into relation (15). Moreover, using expression (14) for the 
partial-wave amplitudes, we arrive at 



A 



£,m 



X 



+ iJ2 f^uPuhf^ 



{PuTq) 



tAv 



V] 



9au 



ti{i 



(19) 
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Equation (19) is the central result of our paper. It explicitly relates the ingredients of elastic 
scattering to the amplitude for production in the spectator approximation. We were able to 
achieve this because in the RSE one can determine in an analytically closed form all terms of 
the perturbation expansions (5) [2] and (15). Hence, relations (11), (12) and (13) can be derived 
and explicitly verified. We may thus conclude that at least for a nonrelativistic (NR) microscopic 
model, i.e., at low energies, production and scattering are related to one another through Eq. (19). 



4.4 P = T/V 

Expression (18) takes an extremely simple form in the case that all inelasticity is either absent 
or neglected. For the £-th partial wave of the production amplitude (15), we then obtain 

cx Xje{pro)Q\:l{E) ^ . (20) 

This is exactly the generic form of the production amplitude used in a paper by Roca, Palomar, 
Oset and Chiang [7], when at the /o(600) resonance the inelastic contribution KK txtx is 
neglected, resulting in a oc T /V . Here, we get in the 1-channel case from Eq. (12) that A = Z, 
which then precisely yields T/V = l/V. 

4.5 The meson-loop phase 

In the one-channel approximation, we obtain from the scattering amplitude (5) for the cotangent 
of the scattering phase shift 5^^\E) the expression 

/ \ rit>{pro) 1 
cotg 5(^)(E) = —— , (21) 

Jiipro) 2X^^ipjj{pro) A^'^ 

where the spherical Neuman function is represented by rti. 

Now, T) in formule (20) is related to A in formula (21) through Eq. (11). After some algebra, 
we get 



«<" oc X„ (pro) (E)\l- /"^fU cos (i">(B)) . (22) 

y J I {pro) I 111 {pro) I ^ ^ 

For 5*- waves {i = 0) this becomes 



With respect to the dependence on the phase S^^\E), this expression has exactly the same form 
as the 5"- wave production amplitude given by Boito and Robilotta in Ref. [41], which is based on 
Watson's formalism [6] via the work of Pennington [42]. For the meson-loop phase u!{s) defined 
in Ref. [41], we obtain here pro. However, our resonance poles are determined in quite a different 
manner than in Ref. [41]. Whereas in the RSE the resonance poles are all contained in A in 
expression (21) for the cotangent of the phase shift, in the formalism employed in Ref. [41] each 
of the resonance poles for S'-wave production has to be put into the corresponding expression by 
hand, one by one. 
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4.6 Breit-Wigner resonances 



Again in the one-channel case, one deduces from Eq. (7) for T) the form 



n=0 



E — E„ 



fJ'Pje {pro) h'-p (pro) 



(24) 



For small X one finds a zero of V in the vicinity of En, say at £"„ + AEn, where 



AEn ^ 



2X' 



F^"^^ (ro) /i„p„ [je (p„ro) rn (p„ro) - i (p„ro) } . (25) 

Here, and are the reduced mass and relative hnear momentum of the two-meson system 
at E = En, respectively. Note that the imaginary part of AEn is negative, as it should be for 
resonance poles in the second Riemann sheet. Below threshold we obtain poles on the real energy 
axis, since ijih^P becomes real for purely imaginary arguments. The latter poles represent two- 
meson bound states, as argued above. For the following discussion we shall only consider poles 
above threshold. 
For V we obtain 

V^'\E) ^X{{E-En-AEn) . (26) 



Consequently, denoting the residue at the n-th pole by a„, we get 

1 an 



oc 



E 



(27) 



pW(^) V {E- En -AEn) 

which is nothing but a Breit-Wigner [43] expansion over a series of resonances, as employed in 
the isobar formalism [40,44-47]. 



4.7 Overlapping resonances 

Of course, things become more involved than in Eq. (27) when X is not small and resonances 
start to overlap. Overlapping resonances have been studied extensively in the past [48]. Here, it 
is no longer possible then to deduce simple approximations for expression (24). 

Besides extending the formalism of Ref. [6] to coupled channels and overlapping resonances, 
our work also seems to interpolate between the results of Ref. [42] and Ref. [7] . 



4.8 The X-matrix 

The i^-matrix, which is related to the tangent (s) of the scattering phase shift (s), is defined by 

K = T [I + zT]"^ . (28) 

As follows from the unitarity condition, K is a, real (symmetric) matrix for real CM energy E. 

In the one-channel approximation and in a particular partial wave, K is given by the inverse 
of expression (21) for the cotangent of the scattering phase shift. For more channels, relations 
like Eq. (21) become very complicated expressions in terms of A and "D. The reason is that the 
inverse of the expression (5) has to be determined. Numerically this is no problem, of course, but 
analytically it is extremely tedious in the general multichannel case. In particular, for a relation 
between the common denominator T) and K, which is needed for the leading term in expression 
(19), nothing simpel follows. Moreover, the pole positions for both scattering and production 
stem from D, and not from K. Hence, the excercise to express the production amplitude in terms 
of the /^-matrix seems pointless. 
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5 Summary and Concluding Remarks 



The two-meson production amplitude (19) has been rigorously calculated, to all orders, from 
a relatively general expression for a two-meson scattering amplitude (Eq. (5)) dominated by s- 
channel resonances. The latter had already been succesfuUy tested for cc and bb states, mesons 
with open charm and bottom, and also in the light-quark sector. 

One might object that a model with no t-channel exchanges is too restricted for drawing 
general conclusions. However, one should be aware of the — quoting Tornqvist [49] — "well- 
known dual-model result for qq resonances that a sum of s-channel resonances also describes 
t- and M-channel phenomena." In the context of duality, Harari [50] formulated a necessary 
condition for an s-channel description to reproduce certain t-channel effects, namely the existence 
of "strong correlations between the different s-channel resonances." Well, this is exactly what 
our infinite RSE sum over confinement states guarantees. Further proof showing the RSE model 
to be realistic is its correct threshold behaviour in elastic vrvr scattering [51]. 

Another possible critique of our method could be its NR nature. Nevertheless, in practical 
phenomenological applications to spectroscopy and elastic scattering, relative momenta and re- 
duced masses in the two-meson channels have been consistently defined in a relativistic way, thus 
ensuring proper kinematics at much higher energies than the underlying NR formalism seems to 
support. Such a minimal treatment of relativity is indeed common practice in many relativised 
quark models. Our successful description of the spectroscopy and scattering properties of the 
light scalar mesons [52] provides additional evidence that this approach to relativity is reasonable. 
This is also supported by our very recent first application of the present production formalism 
in the single-channel case [1]. 

It thus seems fair to conclude that production amplitudes can in general contain terms which 
are not proportional to scattering T-matrix elements and, moreover, that the proportionality 
coefficients are complex. 
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A Generic relation between production and scattering 



In order to arrive at a relation equivalent to Eq. (17) for the here proposed scattering and 
production amplitudes, we define 



ij 



A 



For this object, also using relations (6), one easily finds 



(29) 



A K'-l" A 



(30) 
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Furthermore, we define 



= VfMPijeiPiro) Igai - 2zA^ ^ f^uPu jdP>yro) hf\p^ro) 9au > , (31) 
for whicli, by substituting definition (29), we may also write 

^ai = 9aije{Piro) y/fiiPi + 1 9au ^Jl^-uPu {PuTo) T^^^ • (32) 

For this object we study, in accordance with relation (17), the imaginary part 
= 9au y/fJ'uPu {je [PvTo) 3fte (t/J^) - (p^ro) 3m (t/J^)} , 

(2) (1)* 

where we denote the spherical Hankel function of the second kind hy h\ = h\ = ji — ine. 

Next, we use the fact that 3fJe (T) = T* + i'^m (T), and, moreover, substitute subsequently 
relations (30) and (32): 

(^A';^}) = Y 9au V/VP^ {je {PuTo) T^^ + i h^P {PuTq) (t^^^) } (34) 

= Y 9au y/l^uPu je iPvTo) T^!) + 'i^YY 9av y/fJ'uPu {PuTq) T^^TS 

= Y ^tv ^9auje (Pi.ro) ^/|J,uPu + i Y 9au' y/Jh7p7 (Pu'ro) T^jj 

/ y iv au • 
u 

This demonstrates that for our amplitudes a relation exists which is equivalent to the one shown 
in Eq. (17). 
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